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Energy Dependence Systematics of Strange and Multi-Strange Particle Production 

J. Speltz 1 * for the STAR Collaboration 

* Institut de Recherches Subatomiques, Strasbourg, France 

We present results on systematic measurements of strange and multi-strange particles with 
the STAR detector for center of mass energies per nucleon pair (y^sjvjv) of 62.4 and 200 
GeV in ultra-relativistic Au+Au collisions at RHIC. We use these results to characterize the 
chemical and kinetic freeze-out properties of the fireball produced by the collision. This is 
done by comparison to statistical and hydrodynamical models as well as parameterization. 
We emphasize particularly the energy dependence of these features comparing measurements 
obtained at RHIC with different energies but also at the lower energies available at the SPS. 

PACS: 25.75. -q, 25.75.Dw 

1 Introduction 

The measurement of strangeness in relativistic heavy ion collisions has been suggested to be 
valuable for gaining insight in the created system, as strange valence quarks present in the final 
state of the collision are not existant in the incoming beams. Studying the chemical or kinetic 
freeze-out properties of strangeness in general and compared to non strange particles in particular 
may then reveal different aspects and characteristics of the fireball they originate from. With the 
results from the 62.4 GeV energy that stands between the top RHIC and the SPS energies a more 
complete picture of the excitation function of these properties can be accessed. 

The data used for the presented results are from the STAR experiment for Au+Au collisions 
and were obtained with the main tracking device of STAR [1], the Time Projection Chamber 
(TPC). The 200 GeV strangeness results are from the RUN II data [2], while the 62.4 GeV 
results are from the RUN IV data [3,4]. 

2 Chemical Properties 

The chemical freeze-out marks the end of all inelastic interactions and fixes the relative aboun- 
dance of the different particle species. Statistical models compute these production rates using 
a limited set of parameters and assuming the freeze-out occurs from a statistically equilibrated 
matter. Hence one can consider that these parameters characterize the chemical properties of the 
system. These parameters are the temperature at freeze-out T c / l; the chemical potentials of light 
quarks fi q (q = u, d) and strange quarks fj, s as well as the strange quark phase-space occupancy 
75. For a given collision energy and system, the temperature, potentials and occupancy factor 
are adjusted to the measured "stable" particles ratios (no resonance included). 
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Fig. 1. "Stable" particle ratios as obtained by mea- 
surements for most central events at ^/sjvjv = 
62.4 GeV (circles) together with statistical model 
predictions using [5] (lines). The experimental 
data include statistical and systematic errors. 
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Fig. 2. Centrality dependence of (a) chemical 
freeze-out temperature T c h and (b) strangeness 
saturation factor 7s for 62.4 GeV (closed circles) 
and 200 GeV (open circles). Statistical and sys- 
tematic errors are included. 



FiglO shows the different particle ratios as obtained by measurements for central collisions 
at y/spjN = 62.4 GeV together with a statistical model fit [5] giving T c h = 161 ± 7 MeV and 
Hb = 87 ± 13 MeV. The equivalent figure for 200 GeV can be found in [6]. When comparing 
to the results from SPS [7], it can be assessed that statistical models have been very successful in 
predicting "stable" particle ratios, including multi-strange baryons, on a large range of energies. 

Fig|2] shows the evolution of T c h and 75 with re- 
spect to the mean number of participants (N part ) as 
determined by a Glauber model calculation [8], for 
62.4 and 200 GeV. Firstly, it can be seen that T c h is 
almost constant over the total centrality range and for 
both energies. Its value of approximately 160 MeV is 
close to the predicted LQCD phase transition temper- 
ature of around 170 MeV [9]. Secondly, 75 shows 
a centrality dependence rising from most peripheral 
collisions and reaching unity for most central colli- 
sions. Here again the results from 62.4 GeV and 200 
GeV are very similar and may reveal equilibration of 
strangeness for most central collisions at RHIC ener- 
gies, with contrast to SPS ones where also a central- 
ity dependence for 7s has been manifested, but the 
value for most central events seems to be less than 
unity [10]. 

On Fig[3]the evolution of the strange anti-baryon 
to baryon ratio for central collisions at mid-rapidity with collision energy, from AGS over SPS to 
RHIC, is presented. With the newly added 62.4 GeV value these ratios yield a smooth increase in 
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Fig. 3. Central, mid-rapidity strange anti- 
baryon to baryon ratios as a function of 
ysjvjv [11]. 



Energy Dependence of Strangeness Production At RHIC 



going from a transport preponderant regime at the lower energy to reaching a baryon free regime 
where quark pair production is dominant at the highest RHIC energy. 



3 Dynamical Properties 



While chemical properties can be derived from transverse momentum (pr) integrated values, px 
dependent distributions, such as the particle spectra and differential elliptic flow [12] measure- 
ment are needed to gather information on dynamical characteristics. These are among others 
the kinetic freeze-out temperature Tuin at which all elastic interactions end, but also the flow 
velocity (elliptic and radial) developed during the expansion of the medium. 
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Fig. 4. One and two a contours in Tk in — (fir) 
space obtained by blast-wave parameterization for 
most central collisions at different energies [2, 16] 
and for different particle species. 
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Q~ spectra obtained by measurement at y / ijvjv = 
62.4 GeV (black points) and comparison to blast- 
wave parameterization (black lines) and hydrody- 
namical calculation [13] (grey lines). Errors on the 
spectra are statistical only. See text for more details. 



It has been tried to infer information on these parameters by using either an ideal hydrody- 
namical description [13] which gives a prediction on the spectra and elliptic flow or a hydrody- 
namical inspired blast-wave parameterization [ 14, 15] which permits the extraction of parameters 
by doing a direct fit on the spectra. In the case of the blast-wave these parameters are T^ n , the 
mean tranverse flow velocity (/3t), but also the velocity radial profile n. In the used hydrodynam- 
ical model, in addition to the decoupling temperature Td ec which is the average temperature on 
the decoupling surface, a variable that parameterizes the initial transverse boost at thermalization 
(noted a) was used. 

Fig@] shows confidence-level contours in T^i n — ((3t) space as obtained by blast-wave fits 
on the spectra of different particle species and different energies. The 62.4 GeV and 200 GeV 
contours reveal a clear discrepancy of the transverse flow of the light particles (7r, K and p) and 
the multi-strange baryons (H and f2). Such a disagreement may also be seen on T^n at 200 GeV, 
but is not as obvious at 62.4 GeV. It has been argued at 200 GeV [2] that Tkin is compatible 
with T c h for multi-strange baryons and thus leading to limited interactions of these particles 
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after chemical freeze-out. At 62.4 GeV T c h is comparable to the 200 GeV value whereas Tkin 
obtained by blast-wave fits for H and SI is lower as at 200 GeV. This would then lead to the 
possibility that the multi-strange baryons may continue interacting a little longer after chemical 
freeze-out at the lower energy than at the higher energy. Additionally it can also be seen that (fir) 
is increasing with increasing energy for the light particles, while there may be an indication for 
an increase of temperature for the multi-strange baryons when going from 62.4 GeV to 200 GeV, 
although the values are still compatible within errors. 

Fig|5] shows the fl~ spectra obtained by measurement at 62.4 GeV together with results 
from blast-wave parameterization (black lines) and hydrodynamical calculation (grey lines). 
The black dashed line shows the best fit from blast-wave with the data. This best fit gives 
Tkin ~ 120 ± 30 MeV and ((3t) ~ 0.47 ± 0.05 c. These values show a little dependence 
of the velocity profile which was fixed ton = 1. The black dotted line gives the spectra obtained 
by blast-wave when fixing the parameters to the values obtained by a simultaneous fit to the 
spectra from ir~~, ir + , K~, K + , p and p (Tkin — 90 MeV and ((3t) = 0.57 c). As has already 
been observed qualitatively on Fig[4j the spectral shape of the black dotted line is in clear dis- 
agreement with the measurement. The difference in the lower pr part may partially be attributed 
to incompatibility in flow, while the deviation in the higher pr part may be due to discrepancy 
in freeze-out temperature between the light particles and the Q. As already indicated in Fig|4]the 
former seems to be more pronounced. 

This observation is supported by the spectral shapes obtained by hydrodynamical calculation. 
The use of a temperature Tdec = 100 MeV (grey dotted line on Fig0 and a = 0.02 fm _1 (the 
use of a = fm _1 reveals a larger deviation from measurement at higher pr) gives a form that 
is in agreement with the data, while the shape obtained with a higher temperature (grey dashed 
line) clearly misses that from the data. This has also already been seen at 200 GeV [13]. 

We discuss here only the shape of the spectra and not its normalization as the used hydrody- 
namical model does not include a strange chemical potential that may have significant influence 
on the strange particle multiplicity. Therefore the normalization has been scaled in order to re- 
produce the data. Additionally systematic differences were checked between the hydrodynamical 
and blast-wave temperature by producing a spectra with hydrodynamical model using a certain 
temperature (Td ec ) and then performing a blast-wave fit. We obtain a good agreement, however 
the values obtained by the blast-wave fit are systematically lower than the hydrodynamical tem- 
perature with a smaller difference for lower Td ec (3 — 10 MeV at Td ec ~ 100 MeV) than at higher 
T dec (15 - 25 MeV at T dec » 160 MeV). 

Finally, we also investigated differential elliptic flow resulting from the azimuthal asymmetry 
in momentum space. This elliptic flow has been compared to hydrodynamical prediction. The 
62.4 GeV results are not shown, but give similar results than at 200 GeV [17]. The mass ordering 
at 1ow-pt is well reproduced by hydrodynamics as well as the magnitude up to px of around 
2 GeV/c. This was not the case at the SPS, where the hydrodynamic limit was not reached [18]. 

4 Conclusion 

The high quality and quantity of data accumulated both at SPS and RHIC for different systems 
and at different energies allow for systematic comparison and excitation function of physics ob- 
servables including strange and even multi-strange particles. 62.4 GeV data give qualitatively 
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comparable results as the 200 GeV measurements in term of chemical and dynamical properties. 
The hydrodynamical calculation and blast-wave parameterization seem to work as good at the 
lower RHIC energy (62.4 GeV) as at 200 GeV and give comparable results for freeze-out tem- 
perature. However, the transverse flow obtained by blast-wave for the multi-strange particles is 
significantly lower than that for the light particles. Finally, although ideal hydrodynamics seem 
to work quite well, it is still interesting to see how models that are not assuming such an ideal 
behavior can describe the data, in order to get the entire picture of the properties of the matter 
created at RHIC. 
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